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ABSTRACT 

Aims. We present an accurate characterisation of the high-resolution X-ray spectrum of the Narrow Line Seyfert 1 galaxy Arakelian 
564 and put it in to context with other objects of its type by making a detailed comparison of their spectra. 

Methods. The data are taken from 5 observations with the XMM-Newton Reflection Grating Spectrometer and fitted with various 
spectral models. 

Results. The best fit to the data identifies five significant emission lines at 18.9, 22.1, 24.7, 29.0 and 33.5A due O VIII Lya, O Vll(f), 

N VII Lyo-, N VI(i) and C VI Lyo- respectively. These have an RMS velocity of ~ 1100 kms"' and a flow velocity of 600 kms"', 

except for the O Vll(f) emission line, which has a flow velocity consistent with zero. Two separate emitting regions are identified. 
Three separate phases of photoionized. X-ray absorbing gas are included in the fit with ionization parameters log^ = -0.86, 0.87, 
2.56 and column densities Njj = 0.89, 2.41, 6.03 x 10^" cm"- respectively. All three phases show this to be an unusually low velocity 
outflow (-10 + 100 kms"') for a narrow line Seyfert 1. We present the hypothesis that the BLR is the source of the NLR and warm 
absorber, and examine optical and UV images from the XMM-Newton Optical Monitor to relate our findings to the characteristics of 
the host galaxy. 

Key words. Galaxies: active - Galaxies: Seyfert - Galaxies: individual: Ark 564 - X-rays: galaxies - Techniques: spectroscopic - 
Line: identification 



' 1. Introduction 

Narrow Line Seyfert 1 galaxies (NLSls) were initially classified 
by Osterbrock & Pogge (1985) and are defined as having H/3 
FWHM < 2000km s '. Their X-ray spectra often display more 
rapid variability (e.g. Boiler et al. 1996; Laor et al. 1994) and 
are steeper in the soft X-ray band than those of 'normal' broad 
line Seyfert Is (e.g. Leighly 1999). From the ROSAT All Sky 
. Survey it was found that approximately half of the sources 
' in soft X-ray selected samples of AGN are NLSls (Grupe 
', 1996; Hasinger 1997). They are believed to be 'high-state' 
active galaxies, with low black hole masses for their luminosity, 
and so with high accretion rates relative to Eddington (Pounds 
et al. 1995; Boroson 2002). 

Observations indicate that the majority of Seyfert Is dis- 
play evidence for 'warm absorbers' (Blustin et al. 2005), i.e. 
clouds of photoionized X-ray absorbing gas in our line of sight 
(e.g. Komossa 2000), yet still very little is known of the ab- 
sorbers origin, location or structure. The warm absorber is often 
characterised by its ionization parameter, ^. This is a measure of 
the ionization state of the material and is defined as 

L 

where L is the source luminosity (in erg s"' ), n is the gas number 
density (in cm"^) and r is the distance of the absorber from the 
central engine (in cm) (Tarter et al. 1969). 
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Arakelian 564 is the brightest known NLSl in the 2-lOkeV 
range (L2-Wkev — 

2.4 X 10^3 ergs-'. Turner et al. 2001). 
This, coupled with its close proximity, z=0. 02468 (Huchra 
et al. 1999), and relatively low Galactic column density, 
Nh = 5.34 X 10^"cm-2 (Kalberla et al. 2005), make it a very 
interesting target to investigate. This object has been studied 
across all wavebands (e.g. Shemmer et al. 2001; Romano et al. 
2004). Here we have merged data from 5 observations with the 
XMM-Newton Reflection Grating Spectrometer (RGS), allowing 
us to explore, in unprecedented detail, the high-resolution X-ray 
spectrum of this object. 

iMatsumoto et al.1 (|2004 report on a 50 ks Chandra HETGS 
observation. They fit the hard X-ray spectrum with a power law 
of photon index of 2.56 + 0.06 and add a blackbody, of temper- 
ature 0.124 + 0.003 keV, to fit the soft excess. They detect an 
edge-like feature at 0.712keV, and two phases of absorption by 
photoionized gas (phase 1: log ^ ~ 1, Nh ~ 10^' cm"^; phase 2: 
log^ ~ 2, Nh ~ 10^' cm"2). Vig naU et all (2004) have analysed 
two XMM-Newton observations from June 2000 and June 2001. 
They fit the spectrum with a steep power law (F - 2.50 - 2.55) 
plus a soft blackbody component (kT~140- 150 eV). They also 
identify an edge-like feature in the EPIC data at ~ 0.73 keV, 
which they interpret as an OVII K absorption edge, and find 
evidenc e for a broad iron emi ssion line, at ~ 1 keV, in one obser- 
vation. I Crenshaw et al.l (1200 2) have analysed the UV spectrum 
using data from the STIS onboard the Hubble Space Telescope. 
They identify UV absorption lines due to Lyo-, N V, C IV, Si IV 
and Si III centred at a radial velocity of -190kms ' relative to 
the systemic velocity. They identify this with a dusty lukewarm 
absorber with column density Nh = 1.62 x 10^' cm"- and 
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ionization U=0.033 (this ionization is equivalent to log^ ^ 1.2, 
converted using the method of George et al. 1998). 



iPapadakis et al.l (l2007h analysed the most recent (2005), and 
longest (~ 100 ks), observation with XMM-Newton, focusing 
on the EPIC data. They find that the soft excess cannot be 
parametrized either by a multiple power law or by a power law 
plus a single blackbody, in contrast to previous findings. They 
fit the MOS and FN spectra with either a power law of slope 
2.43 + 0.03 and two blackbodies, with kT~0.15 and -0.07 keV, 
or with a relativistically blurred photoionized disk reflection 
model. They also find evidence for two phases of photoionized 
X-ray absorbing gas with ionization parameters log^ ~ 1 and 
~ 2 and column densities N/^ ~ 2 and 5 x 10^" cm"^. 



Th e same 2005 data were also analysed by iDewangan et al.l 
(|2007|) . They concentrate their analysis on the EPIC data, but 
also study the RGS spectrum and find two warm absorber 
phases with ionization parameters log^ ~ 2 and < 0.3, column 
densities Nh ~ 4 x 10^" and 2 x lO^^cm"^ and flow velocities 
V 300 and -1000 kms"' respectively. 

Here we investigate the combined RGS spectrum of 
Arakelian 564. Section|2]gives information on the individual ob- 
servations and describes our data reduction process. The spectral 
fitting and best fit results are presented in Section [3] A full dis- 
cussion of these results and comparison with other observations 
is given in SectionH] 



2. Observations and Data Reduction 

The majority of our data were obtained with the Reflection 
Grating Spectrometer (den Herder et al. 2001) onboard the 
XMM-Newton observatory. This spectrometer provides un- 
paralleled sensitivity in the 6 -38 A (0.3-2.5keV) range 
with an approximately constant wavelength resolution of 
70 mA FWHM. The data presented comprise four short expo- 
sure (18, 34, 10, 16 ks) observations, two in June 2000 and two 
in June 2001 respectively, and one long exposure (~100ks) 
observation in January 2005, t h e EPI C data of which have been 
presented by IPapadakis et alj (l2007h . The overall appearance 
of the continuum shape is not significantly different between 
observations; as can be seen in Fig.[T]there is a flux variation of 
around 50% between observations. The X-ray flux is known to 
vary by >50% within individual observations (e.g. Papadakis 
et al. 2007; Vaughan et al. 1999), hence the variations between 
our observations are relatively unimportant. There are also 
no significant changes in the warm absorber features between 
observations as we shall demonstrate in Section U4l Therefore, 
data from all five XMM-Newton observations can be combined 
to give the best possible signal-to-noise spectrum. 
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Fig. 1. All five XMM-Newton RGS observations of Arakelian 
564 plotted separately. A similar spectral shape is seen for all 
observations. The two high luminosity spectra are from the 2000 
observations, the three lower are the two 2001 and one 2005 ob- 
servations. 
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Fig. 3. A graph of Ax^ against wavelength to show the signif- 
icant spectral lines. A A^^ of 9 implies a 3cr significance. The 
positive values correspond to emission features and the negative 
to absorption features. The strongest absorption lines are identi- 
fied. 



for use with high resolution X-ray spectra. The data presented 
in the figures have not been corrected for redshift or Galactic 
absorption unless specified otherwise. 



The data were processed using the XMM-Newton Science 
Analysis System (SAS) version 7.0. First and second order 
spectra from both RGS, and all observations, were extracted. 
The spectra and response matrices were resampled and coadded 
to produce a single spec trum (see F ig. |2]) and a single response 
matrix as described in iPage et all (^03). To improve signal 
to noise the data were grouped by a factor of 3, resulting in 
a spectrum with 1000 channels ~30mA wide, well sampled 
with respect to the RGS resolution of ~70mA FWHM. The 
main software used to model the data was spex 2.00, a spectral 
fitting package created by J. S. Kaastra, which is optimised 



In addition to the RGS spectra we also investigate the na- 
ture of the host galaxy using archive B-band and I-band images, 
from 1.2ks exposures in November 1998 with the Im Nickel 
Telescope at the Lick Observatory (Schmitt & Kinney 2000), and 
combined UV images, from the XMM-Newton Optical Monitor 
(XMM-OM, Mason et al. 2001) taken simultaneously with the 
RGS observations mentioned above. The pipeline processed 
XMM-OM UV sky-images in each filter were aligned and then 
coadded to produce the image of Arakelian 564 for each of the 
UVWl (220-400 nm), UWM2 (200-280 nm) and UVW2 (180- 
260 nm) filters. 
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Fig. 2. The combined Arakelian 564 data from XMM-Newton RGS fitted with a power law, blackbody. Galactic absorption, 3 phases 
of photoionized absorbing gas and five narrow emission hnes due to O VIII Lyo', O Vll(f), N VII Lya, N VI(i) and C VI Lya. The 
bottom panel shows the residuals for this fit. 



3. Data Analysis 

3.1. Continuum and Emission Features 

In order to correctly identify absorption and emission features 
a line-finding program (Page et al. 2003) is used to detect 
significant absorption and emission lines in the spectrum. The 
program uses a sliding cell (~ lA wide) detection technique, 
rejecting all points more than 2cr deviant to create a smoothed 
continuum. A grid search is then performed, adding a Gaussian 
to the continuum model at each point and calculating the 
goodness of fit by minimising (Lampton et al. 1976). This 
produces a graph of against wavelength as shown in Fig. |3] 
A Ax^ value of 9 corresponds to a 3cr significance and so, with 
the Arakelian 564 data, only one spurious statistical fluctuation 
would be expected at this significance level. 

To begin the fitting proc ess our continuum par ameters are 
fixed at those required by the iPapadakis et afl (|2007|) power law 
(f = 2.43) and two blackbody (kT=150, 70 eV) fit. We also 
include Galactic absorption (Nh - 5.34 x 10^° cm"^, Kalberla 
et al. 2005). 

Using Fig. [3] as an initial guide, we can fit a number of 
emission lines to the data. The strongest three emission lines at 
~ 19.40, 25.35 and 29.75A are well modelled with Gaussian 
profiles, and match the wavelengths of O VIII Lya, N VII Lyo- 
and N VI(i) respectively. Using the Gaussian component in spex 
these emission lines are added to the above continuum and the 
profiles fitted to the data. The fit can be seen in Fig.|4]on an ex- 



panded scale. Also included in the fit are emission lines from 
OVII(f) and CVI Lya at ~ 22.65 and 34.35A respectively, 
which are required by the data (improvement in^^ of ~ 50 for 6 
degrees of freedom: normalisation, wavelength and FWHM for 
each line) but are not as significant in Fig. [3] due to the close 
proximity of absorption features. The parameters for all emis- 
sion lines can be seen in Table[T] which lists the identification of 
each emission line, FWHM, flux and velocity characteristics. 



3.2. Absorption Features 

The most significant absorption lines in Fig. [3] match the 
wavelengths of NelX, OVUI, Galactic OI and Oil, N VII and 
CVI. There are a number of spectral components in spex that 
can now be included in the fit to characterise these absorption 
lines. Here we use the xabs component. This is a self-consistent 
model, which contains all of the ions present in a photoionized 
gas at a particular ionization level. Free parameters in the fit 
are the ionization parameter, column density, RMS velocity and 
flow velocity; the abundances have been set at the Solar values 
(Verner et al. 1996). Multiple 'phases' of gas can be added to 
the model, using multiple xabs components, to reproduce the 
observed spectrum. We us ed the Arakeli a n 564 spectral energy 
distribution (SED) from iRomano et all (|2004 (fliis SED is 
consistent within errors with the OM and RGS average fluxes 
measured from our observations) to model the ionizing flux and 
compute the absorption profile in spex. 
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Fig. 4. A zoomed in section of the data to show the most prominent absorption and emission lines. 
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Table 1. Table of parameters for the five narrow emission lines included in the fit (errors are 90% confidence). 



Line 


Observed A 
A 


Rest frame A 
A 


FWHM 
A 


FWHM 

km s"' 


Flow Velocity 

km s"' 


Line Luminosity 

lO^Vhs-' 


O VIII Lya 
OVII(f) 
NVIILyor 
NVI(i) 
C VI Lya 


19 41+0.01 
22 65+° °^ 

■:i4+0.02 

70 74+0.04 
'^-0.05 

34 36+°°2 


18 94+0.01 

22 10+°°^ 
94 7^+0.02 

29.02:0;°^ 
a-j c')+0.02 


"•^^-0.22 


1400+^*"° 

1500!?! 
850!™ 
850!^» 


-440+200 
QO+iooo 

'"-700 

-650!2oo 
-630!«« 
-1800!2» 


9 1 1 +0.98 
■^•J-i_0.78 

-'•'^"-1.84 

1 7C+O.47 
'^•'-'-0.70 

0.88!«i^ 
2.26!«Ji 



Table 2. Table of parameters for the three phases of photoionized, X-ray absorbing gas (errors are 90% confidence). 



Phase 


Flow Velocity 


RMS Velocity 


Log^ 


Column Density 




(kms"') 


(kms"') 


(erg cms"') 


(102''cm-2) 


Phase 1 


9/1+200 

^^-'-100 


100+350 


-0.86!«:|« 


89+°-^'* 


Phase 2 


_40+50 


60+2° 

""-10 


0.87!jJ:?^ 


9 41+0.34 
^•^^-0.34 


Phase 3 


_ 10+40 

■^"-100 


100+^° 


^•-'"-0.04 


6 03+''''' 
"•"■-^-0.86 
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The best fit to the data includes three phases of photoionized 
gas, the model profile of which can be seen in Fig. |5j the fitted 
spectrum is in Fig.|2l The parameters of the model are listed in 
Table|2l each phase has a different ionization parameter and col- 
umn density. All three phases have compatible velocity profiles; 
however, the flow velocities are not compatible with those of the 
emission lines (see Table [TJ, which are faster outflowing, with 
the exception of the O Vll(f) line. The uncertainty in the wave- 
length scale of RGS data is ~ 8 mA, giving rise to a system- 
atic uncertainty on the outflow velocities of at most 100 km s"^ 
Therefore, it is unlikely that the absorber and the majority of the 
emission lines are connected. However, Phases 2 and 3 are com- 
patible with the O Vll(f) emission line in ionization and velocity, 
so it is possible that these are connected. 

3.3. Ionization structure 

It is possible that a continuous distribution of gas with a range 
of column densities and ionizations, rather than three discrete 
phases, is what is observed in Arakelian 564. To investigate this 
possibility further the data were fitted with the same continuum 
and emission lines as above, but using the warm component in 
place of the three xabs components, warm uses the equivalent 
of 19 XABS phases by default, each having an incremental 
logarithmic increase in ionization ranging between two user set 
ionization parameters. The column density is a free parameter 
for each of the 19 phases, but all have the same velocity profile. 

The result is similar to that using xabs. The best fit gives 
an RMS velocity value of 70^^kms"', and flow velocity of 
-40^35 km s"', which are compatible with the results from the 
xabs fit. Fig. |6] shows the column density plotted against the 
ionization parameter for the warm component. Three peaks are 
clear in the data, at approximately the same ionization values 
and column densities as the three xabs phases fitted previously. 
This adds to the credibility of the initial fit and suggests that 
there are indeed three discrete phases of photoionized gas, 
rather than a continuous distribution of ionization parameter. 

To investigate the possibility of all three phases co-existing 
in the same location we use the pressure form of the ionization 
parameter, S. This is related to f by 




1 

Ionization {log ^} 



Fig. 6. The distribution of column density versus ionization pa- 
rameter using the warm component in spex to fit the data with a 
continuous phase of gas. 
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Fig. 7. Thermal stability graph showing the equili brium temper- 
ature a s a function of S for the SED described in iRomano et al.l 
(2004). The positions relating to the three warm absorber phases, 
along with their errors, are indicated. 



0.961 X 10^^ 



4-ncr^P 



(2) 



where Lj„„ is the 1-lOOORyd ionizing luminosity, r is the 
distance of the absorber from the ionizing source, P is the 
pressure and T the temperature (Krolik et al. 1981). 



The SED from 'Roma no et al. I (12004 was used to create 
a graph of thermal stability of H against temperature, that is 
shown in Fig. |2l The best fit ionization parameters and their 
errors, from Table|2] are used to mark the regions corresponding 
to the three phases on the graph. In order for the phases to 
be in pressure equilibrium, and so co-exist, they must have 
overlapping values of E. In Fig. |7]it can be seen that the three 
phases do not overlap in S and so cannot co-exist in pressure 
equilibrium. The majority of AGN di splay an uns t able r egion 
where d{lnT) I d{lnE) < 0, as noted in lKrolik etall (Il98ll) . The 
lack of any unstable region is unusual, but the shape of the 
cu rve is very similar to th at seen in the NLS 1 IZw 1 as presented 
bv lCostantini et alj ( l2007h . 



Constraints can be put on the location of the warm absorber 
phases. We can estimate the maximum distance that the gas can 
be from the continuum source, and still be ionized at the mea- 
sured levels. We can assume that the thickness of any gas phase, 
Ar, has to be less than or equal to the distance from the source, 
R: 



Ar 



< 1. 



(3) 



The column density, Nh, as a function of density, n{R), is given 
by 



Nh ~ n(R)C,Ar 



(4) 



where Cy is the volume filling factor. This, together with 
Equation[T]gives the expression for the maximum distance 



R < 



LionCy(R) 

4Nh 



(5) 




(Blustin et. al 2005). We can estim ate L,o„ by integrating over 
the SED from iRomano et alJ (|2004 : LogL,o„ = 45.2. We use 
a typical value of Cy{R) of 0.03 (from the average Cy(R) of 
NLSls calculated by Blustin et al. 2005), since Cy(R) cannot 
be calculated for Arakelian 564 because of its small outflow 
velocity. Values for the ionization parameter and column density 
are taken from Table |2l This gives values of R < lOOOkpc, 
lOkpc and lOOpc for Phases 1, 2 and 3 respectively. 

Additional information can be gained from l ooking 
at t he UV spec t ra pre sented in ICrenshaw et alj (|2002|) 
and iRomano et alj (|2002|) . It is clear from the UV absorp- 
tion and emission line profiles that the absorption eats away 
the emission lines and so the absorption occurs between us and 
the emission line region. Therefore, we can use the distance 
of the UV emission line region as a minimum distance to our 
warm absorbe r for Phase 1 and 2, since these are seen in the 
UV spectrum. iRomano et al.l (l2004i) calculated the size of the 
H/3 emitting line region, using the FWHM of UV and optical 
emission lines, as Rblr - 0.008 + 0.002 pc. 

ICrenshaw et alj (|2002|) find that the UV absorber must lie 
outside the NLR and use the Lyo- spatial profile to put a limit on 
the size of the NLR of < 95 pc. Since Phases 1 and 2 have been 
shown to contribute significantly to the UV spectrum they must 
also lie outside this region. 

We can also estimate the minimum distance to the absorber 
by assuming that the outflow velocity has to be greater than or 
equal to the escape velocity of the black hole. Arakelian 564 
has very low velocity outflows, or the absorber could even be 
inflowing; using the maximum outflow velocity for Phase 3 of 



1 10 km s we can calculate a minimum distance to the absorber 
of 



where G is the gravitational constant and M is the mass of the 
black hole (= 3.7 x 10^ Mq, Matsumoto et al. 2004). This gives 
a value of R > 25 pc. This is consistent with Equation |5] which 
gave a value of < lOOpc for Phase 3. Assuming that this Phase 
must be travelling faster than the escape velocity and R < 100 pc, 
we have a minimum outflow velocity of ~ 60 km s"' . 

3.4. Analysis of Individual Observations 

In order to justify coadding all five spectra we must be con- 
tent that there are no significant spectral changes between the 
observations. Fig. [1] shows that changes in the continuum and 
overall shape are minimal. Due to the poor signal-to-noise of 
many of the spectra it is not possible to fit them all separately. 
However, Fig. [1] indicates two observations at a higher level of 
luminosity, and three at a lower level. Any spectral changes that 
are present would be expected between these two states and so 
the two high-state observations were added together and so were 
the three low-states to improve the signal-to-noise. These two 
spectra were then fitted with the above model, including five 
emission lines and three warm absorberphases. In the high-state, 
three warm absorber phases are found with ionization parame- 
ters logf = -0.3 + 0.5, 1.0 + 0.4,2.2 + 0.4 and column den- 
sities Nh = 0.9 + 0.7, 1.4 + 1.1,7.7 + 2.0 X 10-" cm^^. The 
low-state warm absorbers have ionization parameters log^ = 
-0.5 + 0.4,1.1 + 0.6,2.2 + 2.1 and column densities Nh - 
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Fig. 8. Model graph of the expected UV spectrum, due to Phase 
1 and Phase 2 absorption (Phase 1 in black. Phase 2 in red). 

1.0 + 0.9,2.1 + 1.3,7.2+ 1.7 x lO^^cm-^ Comparing these pa- 
rameters with those in Table |2] we see that all are compatible. 
Therefore there are no significant spectral changes between ob- 
servations and they can be coadded to produce one high signal- 
to-noise spectrum. 

4. Discussion 

4.1. X-ray and UV Absorption 

By combining all XMM-Newton RGS observations of Arakelian 
564 we have identified three phases of photoionized. X-ray 
absorbing gas in its spectrum. The model used for the three 
phases can be seen in Fig. |5] Phase 1 (black) contributes 
absorption from very low ionization iron (Fe II-Fe IV); Phase 2 
(red) contributes the mid-ionization iron M-shell UTA features 
along with those from ions such as NVI, OVI and CVI; the 
high ionization iron, oxygen and neon features are produced 
by Phase 3 (blue). The low ionization of Phases 1 and 2 means 
that that they will contribute significantly to the UV spectrum: 
this can be seen in Fig. |8] where strong carbon and oxygen 
absorption features are visible overposed on the extrapolation of 
the RGS continuum by the xabs model component in spex. 
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MCG-6-30-15 (Sako et al. 2003). On the other hand, what we 
see here is a relatively weak absorber with a very low outflow 
velocity. The low ionization phases at low velocity could be 
explained as interstellar gas. We have shown that Phases 1 and 2 
could be up to lOOOkpc and lOkpc respectively from the core; 
however, it is likely that these phases are somewhat closer in, 
as a substantial column of ionized gas would not be expected 
outside of the body of the host galaxy, and so we can put a limit 
on the distance to the absorber of the size of the galactic bulge. 
From examining the B-band image in Fig. |9]it is clear that the 
visible extent of the bulge is at most 3kpc. At this distance the 
gas kinematics will be dominated by the gravitational potential 
of the host galaxy rather than that of the black hole. Phase 3 is 
highly ionized and so must be part of the AGN system, and yet 
it has a very low flow velocity. For Phase 3 to exceed the escape 
velocity it must have an outflow speed in excess of 60 km s"' . If 
it does not, we are looking at either a stationary warm absorber 
or one with a transverse movement with respect to our line of 
sight. 

By fitting the data with a continuous distribution of pho- 
toionized gas, we have shown that three separate phases are 
likely to exist, rather than a continuous outflow with changing 
ionization and column density. We can also say that the three 
phases are unlikely to co-exist in the same region of the AGN 
system, since they are not in pressure equilibrium. If they did 
exist within the same region the pressure difference would force 
the higher pressure region (Phase 1) into the lower pressure re- 
gions (Phases 2 and 3) and what we are observing would not be 
a stable condition. The phases being separate supports the idea 
of the low ionization phases being interstellar gas in the host 
galaxy and Phase 3 being part of the inner AGN system, with 
a far higher level of ionization. We are also able to put reason- 
able constraints on the location of Phase 3. Our calculations in 
Section [33] show the distance from the central source to Phase 
3 to be 25 < 7? < 100 pc. Using the expression for the distance 
to the inner edge of the torus, ~ 1 x L'^^^^ ^ pc (Krolik & Kriss 
2001), we get r,o,.„5 ~ 3.9 pc. Because of their expected contribu- 
tion to the UV spectrum Phases 1 and 2 are likely to lie outside 
the majority of the NLR, putting them well outside the base of 
the dusty torus, eliminating an accretion disk origin; however, as 
we have seen. Phase 3 could be much closer in and could even 
be connected to the NLR itself. 



Our, more detailed, find ings are in reasonable agreement 
with iDewangan et"an (|2007|) . Phase 3 is similar to their higher 
ionization warm absorber phase in ionization and column 
density. Their lower ionization phase lies between our Phase 
1 and 2 in ionization and column den sity. Our Phases 2 and 3 
match reasonably well the results from^Matsumoto et al.l ( 20041) 
in ion ization and velocity profile but the Matsumot o et al.l 
(12004 outflow column densities are larger (phase 1: log^ ~ 1, 



Nh ~ 10^' cm-^; phase 2; log^ ~ 2, Nh 
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21 , 



^). Our 



Phase 2 is al so compatible with the UV absorption previousl y 
identified by ICrenshaw et al.l (|2002|) and iRomano et al. I (I2002h . 
who found a phase with log^ ~ 1.2 and log Nh ~ 10^'cm"^. 
Therefore, it is likely that this and our Phase 2 are the same 
absorber. Phase 3 is too highly ionized to be detectable in the UV. 

The low outflow velocities of these phases of absorbing gas 
are unusual. NLSls are thought to be very luminous for their 
mass, and so in general would be expected to have strong, fast 
outflows such as those seen in NGC 4051 (Ogle et al. 2004), and 



4.2. Emission Features and Connections witli tine Absorber 

It is clear that the X-ray absorption does not originate in the 
same region as the majority of the emission lines. The emission 
lines display higher outflow velocities (see Tables [T] and |2]i, 
suggesting an origin closer to the centre of the system. This 
agrees with the UV data (Crenshaw et al. 2002, Romano et 
al. 2002) in which the absorption lines eat into the emission 
lines, implying that the absorption occurs further out from 
the nucleus than the emission. It is important to note that the 
velocity profiles of the emission lines may be altered due to the 
presence of absorption lines close by, as can be seen in Fig. |4] 
From Fig. |4]it is clear, particularly in the O VIII Lya profile, 
that the absorption lines lie to the red of the emission lines. 
Therefore, the emission line producing gas must be outflowing 
faster than the absorption line region. The O VIII Lya, N VI(i) 
and NVII Lya emission line flow velocities are consistent, 
within errors, indicating a likely common origin. The C VI Lya 
emission line appears to have a far higher outflow velocity than 
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Fig. 9. From left to right: I-band and B-band images, from the Lick Observatory; UVWl, UVM2 and UVW2 images, from the 
XMM-Newton OM, of Arakelian 564. 




Fig. 10. This artist's impression of the core of Arakelian 564 shows the limits on the locations of the broad emission, narrow emission 
and absorption line regions. 



the others. This could be an artifact of absorption feature con- 
tamination, being the least significant of the fitted emission lines. 



O VIII Lya, N VI(i) and N VII Lyo- have velocity widths of 
~ 1000 km s"' making their likely origin the broad line region 
close in to the centre of the system. This would be expected 
to be a region of high density, as is derived for the optical 
and U V broad line region (n^ ~ 10"cm"^) by iRomano et al.l 
(12004 . The presence of an N VI(i) emission line without any 
clear identification of a forbidden line implies a region of 
hi gh density. Follow i ng th e analysis of photoionized plasmas 
bv lPorquet & Dubaul (l2000t) we define the ratio of the forbidden 
line (z) to the intercombination lines (jc + y) as 7? = zlix + y). 
From the data we fi nd a value of << 1 for t he N VI producing 
region . From the iPorquet & Dubaul (l2000l) and iGodet et al] 
(|2004|) plots of the ratio, R, against density, n^, we see that 



R « I for N VI implies a density of >> 3 x 10 
is consistent w ith the BLR density of 
bv lRomano"et al.. (.2004,) . 



iOcm-3. 



This 



The OVII(f) emission line is seen to have a very low 
outflow velocity compared to the other soft X-ray emis- 
sion lines identified here. This could be explained by it be- 
ing produced in a different region of the AGN: the nar- 
row line region. The measured velocities are in agree- 
ment with the optical and UV narrow line region values 
(FWHM< lO CIOkms ', blues h ifted b y 400-1200km s -'), as 
me asured by ICrenshaw et al.1 (|2002|) . iRomano etakl (l2002h 
and lContini et al.l ( l2003l) . Interestingly the profile of the O Vll(f) 
emission fine is also consistent in velocity and ionization with 
Phase 2 or 3 absorption. Therefore, it is possible that a com- 
ponent of the absorber and this emitter have the same origin. 
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The O VII(f) emission line is d ensity-sensitive, so we again 
use the iPorquet & Dubaul ( l200d) analysis. We find a value of 
^ > 3 for the O VII emission triplet. This implies a density of 
He < lO' cm~^, in agreement with NLR conditions and not those 
of the BLR, adding credibility to the theory that the O Vll(f) fine 
originates in a different region to the other identified soft X-ray 
emission lines. 

4.3. Comparison with warm absorbers in other NLS1s 

We have compared these results with a number of other NLSls. 
NGC 4051 displays a similar continuum shape to Arakelian 
564, with a steep soft X-ray continuum fitted with a power 
law of slope r=2.2 (Ogle et al. 2004). These authors note the 
existence of a broad C VI Lya emission line, the core of which 
is found to have a FWHM velocity of 1200 km s"', from a 
region of density = 1 x lO^^cm""' at ~ 2 x lO^'^pc from 
the central source, with a possible outflow in the outer BLR; 
this FWHM is compatible with our 'broad' emission lines, 
though the density is slightly lower This line also has a broad 
component implying far higher FWHM velocities and densities, 
which is not seen in the emission l ine pr ofiles of Arakelian 
564. Also identified, by ' Ogle et all (|2004 . are a number of 
narrow emission lines, due to NVI, OVII, NelX and Si XIII 
with densities n,, = 10^ - 10^ cm"-', consistent with the density 
we find for the narrow line region in Arakelian 564. There is a 
large range in ionization of the absorbers of NGC 405 1 which 
is inc onsistent w ith being in pressure equilibrium, as also seen 
here. lOgle et all (^004) indicate that the narrow lines and an 
ionized absorber in NGC 405 1 may arise in the same region. 

The NLSls MCG-6-30-15 and Markarian 766 also show 
strong emission and narrow absorption lines in their RGS soft 
X-ray spectra (Branduardi-Raymont et al. 2001). Their fits 
include steep power laws (r=2.14 and 2.53 respectively) and 
relativistic emission lines from O VIII, N VII and C VI, as well 
as narrow absorp tion Unes, i n orde r to explain the peculiar shape 
of the spectrum. Lee et al.] ( l2003l) have also analysed Chandra 
HETGS data of MCG-6-30-15 and favour a fit of the spectrum 
using purely a dusty warm absorber model. The analysis of the 
object is still the su bject of much con troversv; ISako et al.l (l2003l) 
have shown that the iLee et al] (l2003h model is a poor fit outside 
of the HETG range and have also made a detail ed analysis 
of the warm absorber present in the two sources. ISakp et al.J 
(l2003l) identify two warm absorber phases in MCG-6-30-15 
with column densities of ~ 2 x 10^' cm"^ and outfl ow velocities 
of ~ -150 and ~ -1900kms ' respectively. iMason et al.l 
(12003') use a longer RGS exposure to fit the absorption lines 
of Markarian 766 individually with column densities of 
3.4 x lO'^-cm ^ for CVI and 2.2 x lO^'^cm-^ for OVII with 
a velocity width of ~ lOOkms ' and a small outflow velocity 
< 100 km S-' 

Two warm absorber phases are also found in the spectrum of 
IZwl by lCostantini et al](l2007h . The two absorbers, which have 
ionization parameters of log ^ ~ and log^ ~ 2.5, have similar 
column densities of ~ 1.3 x 10^' cm"^. The lower ionization 
phase has an outflow velocity of ~ 1700kms"', and the higher 
of < 800 km s"'; implying far higher outflow velocities than 
those in the Arakelian 564 absorbers. The warm absorber phases 
in IZwl are similar in ionization to our Phases 1 and 3; the 
lower ionization phase of the two is contributing significantly 
to the UV spectrum, as seen in Arakelian 564. The shape of the 
thermal stability graph in IZwl is also very similar to that seen 



here (Fig.|7]l in that there is no unstable region. ICostantini et al.l 
(2007) suggest that this could be due to the steep soft X-ray 
spectra displayed by many NLSls. They also find that, like ours, 
their warm absorber phases are not in pressure equilibrium and 
so would require, for example, magnetic confinement in order 
to co-exist. 

Markarian 359 is a N LSl which, lik e Arakelian 564, has a 
very weak X-ray absorber. lO'Brien et alj (1200 lb detect evidence 
for a weak Iron M-shell UTA, but no other significant absorption 
features. Also identified in the RGS spectrum, are three narrow 
emission lines at 13.6, 22.0 and 18. 9A, corresponding to the 
NelX and OVII triplets and the O VIII Lya emission line. The 
signal-to-noise of the spectrum is such that accurate velocity 
values cannot be well constrained. Ton SI 80 has a weak warm 
absorber, thoug h with a larger c olumn than that identified in 
Arakelian 564. iRozariska et alj ( |2004|) identify a number of 
narrow absorption lines with a minimum column density of 
5.1 X 10^' cm"^. There is also emission from CVI, OVIII and 
NeX. Unfortunately, due to the weakness of these features, 
very little information can be determined in order to make a 
comparison with Arakelian 564. 

By comparing Arakelian 564 with other NLSls we can see 
that it is an unusual object, but at least two other galaxies, 
Markarian 359 and Ton SI 80, may have similar warm absorber 
properties. The warm absorber in Arakelian 564 is very weak 
and has a very low outflow velocity. Unfortunately the available 
data for Markarian 359 and Ton S 1 80 are not of sufficient quality 
to determine outflow velocities; the only NLSl warm absorber 
that we have found in the literature, with well defined parameters 
and such low velocities as Arakelian 564, is that in Markarian 
766 (e.g. Mason et al. 2003). It is also unusual to find strong, 
broad, blueshifted emission lines in the soft X-ray spectrum, 
with similar profiles as seen in the UV lines, raising questions 
about what happens to this gas and how it may be connected to 
the narrow emission and absorption line regions. 

4.4. Physical Structure 

We can see from the images in Fig. |9] that the host galaxy is of 
SBb type (e.g. Petrov et al. 1999). The UV images show that 
the galaxy has some unusual characteristics, such as its square 
shape. Where the I-band and B-band images show evidence of 
a fairly typical SBb galaxy, the UV images show some features 
that do not conform to the spiral structure of the galaxy. These 
features are most noticeable in the north-east and south-east 
corners of the UVWl image as indicated in Fig. [TT] They are 
a possible indication of some disturbance to the galaxy in the 
past, producing tidal streams. It is possible that the disturbed 
nature of the host galaxy would be the cause of some of the 
absorption and extinction features observed in Arakelian 564. 

We can estimate the viewing angle of the host galaxy by 
measuring its ellipticity. This is done by measuring the major 
and minor axis of a number of isophotes fitted to the I-band 
image of the galaxy. From averaging these values we find an 
ellipticity value, defined as e = 1 -b/a (a = major axis, b=minor 
axis), of e = 0. 184; this compares well wi th the measured value 
of 0.198 from 'Schmit t & KinnevI ([2000). Using our value of 
e - 0.184, and assuming an intrinsically circular galaxy disk, 
we calculate an inclination of 35 + 10° (errors calculated at 
90% confidence from measuring the ellipticity of the galaxy at a 
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Fig. 11. UVWl image of Arakelian 564 with the contrast set such that disturbance features are enhanced. 



number of different isophotes) in agree r nent w ith the 37° value 
calculated from the ISchmitt & Kinnevl (|2000|) ellipticity. If the 
plane of the AGN accretion disk was aligned with the plane of 
the host galaxy then we would expect to find that the position 
angle (P.A.) of the radio jet in this objec t lies alo n g the same 
axis as the host galaxy minor axis. Sch mitt et al.l (1200 lb find 
that the radio structure of Arakelian 564 lies in a north-south 
direction (P.A. ~ 6°). From Fig.|9]we can see that the P.A. of the 
minor axis of the host galaxy is ~ 50°, and so the two are not 
aligned. 

From our analysis and previous studies of the UV and 
optical spectra it is clear that the BLR in Arakelian 564 has a 
reasonably large outflow velocity of ~ 1000 km s"', the NLR 
a slightly lower one at a few hundred kms ' and the warm 
absorber is outflowing at a maximum speed of ~ 100 kms '. 
All thr ee com ponents d isplay a very similar range of ion- 
ization. ICren shaw et al\ (|2002|) state that the warm absorber 
contains a large amount of dust iE(B - V) - 0.14); using 
the standard Galactic gas-to-dust ratio this implies gas with a 
column density > 8 x 10^" cm"^. At this gas-to-dust ratio, all 
three of our absorber phases would be needed to contain this 
dust and therefore the majority of this material must lie outside 
of the NLR in order to significantly redden the emission lines. 
Warm absorbers, in general, do not contain large amounts of 
dust, making this object quite unusual. 

We can calculate mass outflow rates for the different regions 
of the AGN. For the BLR we first need to estimate the filling 
factor using the equation 



Filling factor ■ 



3L(Ha) 



(7) 



from iDevereux & Shearerl (l2007h . where L(Ha) - 
5.8 X 10^° erg s^' is the luminosit y of the H a line , 
«H — 10^' cm~^ is the density from iRomano et all (|2004|) . 



Ho 



= 9 X 10 



-14 -3 -1 

cm ^ s 



is the effective recombination 



coefficient as defined in lOsterbrock & Ferlandl (12006) . and r is 
the distance of the BLR from the centre of the system. 

This gives a filling factor of approximately 3 x 10"^. 
Using this we can estimate a mass outflow rate o f the BLR as 
~ 1.3 X 10 "^Moyr Foflowing the method of iBlustin et alj 



(l2005h . we calculate the total mass outflow rate of the absorbing 
regions as ~ 0.6, 0.2 and O.lMoyr"' for Phases 1, 2 and 3 
respectively. This is significantly larger than the BLR mass 
outflow rate, which could be explained by the BLR sweeping 
material off the torus and creating the dusty warm absorber. 
Mass loading of the gas would also slow it down, in order to 
conserve momentum, to the required speeds observed in the 
absorber 

Combining this evidence suggests that the BLR, NLR and 
absorber may have a common origin. The outflowing BLR 
becomes more diffuse with greater distance from the central 
system, keeping the ionization of the gas approximately the 
same. Part of this outflow sweeps material off the edge of 
the torus, slowing the outflow and increasing its dust content. 
This then becomes the less clumpy NLR and the diffuse warm 
absorber, as is shown in Fig.fTOl 

Alternatively one or more of the phases of ionized absorber 
could be up to a few kpc from the core. This could be ex- 
plained by the observation that the majority of NLSls contain 
bars (Crenshaw et al. 2003), which enable the transport of dust 
and gas along the leading edge of the bar (Sparke & Gallagher 
2000). This material infalls towards the centre until it reaches a 
stable circular orbit around the galactic bulge. This then either 
stays in orbit or infalls in a circular or barlike motion. This ring 
of gas and dust could be the source of our low ionization ab- 
sorber phases, explaining low flow velocity. Also, disturbance 
to this galaxy, of which we have tentative evidence, could have 
caused gas and dust to be repositioned throughout the system, 
giving rise to dust and gas lanes within our line of sight. 



5. Conclusion 

We present a detailed analysis of the photoionized. X-ray ab- 
sorbing and emitting gas present in the NLS 1 galaxy Arakelian 
564. The high-resolution X-ray spectra are from five observa- 
tions with XMM-Newton's Reflection Grating Spectrometer. 

The absorption profile is fitted with three phases of X- 
ray absorbing gas with average RMS velocity ~ 60 kms"' 
and flow velocity ~ -lOkms Each phase has a different 
ionization (log ^ = -0.86,0.87,2.56) and column density 
(Nh = 0.89 X 10'2",2.41 X 102",6.03 x lO^^'cm^^ respectively). 
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Note that one continuous distribution of gas is unlikely, but not 
ruled out. The low flow velocity of the warm absorber phases 
can be partially explained by assuming the low ionization phases 
to be interstellar gas. The lower ionization phases are expected 
to contribute significantly to the UV spectrum, implying that the 
UV and X-ray absorbers are connected. 

The best fit to the RGS data requires five significant emis- 
sion lines due to O Vlll Lya (18.9A), O Vll(f) (22.1A), NVI(i) 
(29.0A), N VII Lya (24.7A) and C VI Lya (33.5A). By examin- 
ing the outflow velocities and density restrictions on these emis- 
sion lines we have shown that they originate in two separate re- 
gions, the BLR and the NLR. It is possible that these regions and 
the absorber are connected. The outflowing broad line region ex- 
pands as it travels out from the core, maintaining a similar ion- 
ization level. This gas then picks up dust from the torus, and 
slows down. As the gas continues to expand it then displays the 
characteristics of the NLR and the diffuse, dusty warm absorber. 
UV observations suggest a disturbed host galaxy. A tidal stream 
or dust and gas accumulating in the iimer regions of the host 
galaxy offer an alternative origin for the low ionization absorber. 
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